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ABSTRACT 

We present stellar velocity dispersion (cr*) measurements for a significant sample of 40 broad-line (Type 1) 
active galaxies for use in testing the well-known relation between black hole mass and stellar velocity dispersion. 
The objects are selected to contain Ca II triplet. Mg lb triplet, and Ca Hh-K stellar absorption features in their 
optical spectra so that we may use them to perform extensive tests of the systematic biases introduced by both 
template mismatch and contamination from the active galactic nucleus (AGN). We use the Ca II triplet as a 
benchmark to evaluate the utility of the other spectral regions in the presence of AGN contamination. Broad Ee II 
emission, extending from ~ 5050-5520 A, in combination with narrow coronal emission lines, can seriously bias 
cr* measurements from the Mg lb region, highlighting the need for extreme caution in its use. However, we argue 
that at luminosities constituting a moderate fraction of the Eddington limit, when the Ee II lines are both weak 
and smooth relative to the stellar lines, it is possible to derive meaningful measurements with careful selection 
of the fitting region. In particular, to avoid the contamination of coronal lines, we advocate the use of the region 
5250-5820 A, which is rich in Ee absorption features. At higher AGN contaminations, the Ca Hh-K region may 
provide the only recourse for estimating cr*. These features are notoriously unreliable, due to a strong dependence 
on spectral type, a steep local continuum, and large intrinsic broadening. Indeed, we find a strong systematic 
trend in comparisons of Ca Hh-K with other spectral regions. Luckily the offset is well-described by a simple 
linear fit as a function of cr*, which enables us to remove the bias, and thus extract unbiased cr* measurements 
from this region. We lay the groundwork for an extensive comparison between black hole mass and bulge velocity 
dispersion in active galaxies, as described in a companion paper by Greene & Ho. 

Subject headings: galaxies: active — galaxies: kinematics and dynamics — galaxies: nuclei — galaxies: Seyfert 


1. INTRODUCTION 

It has long been known that some galaxies harbor super- 
massive black holes (BHs) at their centers, whose accretion- 
powered luminosity may outshine their entire host galaxy 
(e.g., Lynden-Bell 1969). More recently, we have learned that 
most, if not all, galaxies with a bulge contain central BHs (Ko- 
rmendy & Richstone 1995; Magorrian et al. 1998; Ho 1999), 
whose masses are tightly correlated with the stellar velocity dis¬ 
persion (ct*) of the bulge (the Mbh - cr* relation: Eerrarese & 
Merritt 2000; Gebhardt et al. 2000a; Tremaine et al. 2002). It is 
likely that the Mbh “ cr» relation is established during the active 
galactic nucleus (AGN) phase of a galaxy’s life-cycle, since en¬ 
ergy emitted by the BH may simultaneously limit the gas sup¬ 
ply for building both the bulge and the BH itself (e.g.. Silk & 
Rees 1998). If we are to understand this critical stage in galaxy 
evolution, we must have robust methods to estimate both BH 
masses and cr* in active galaxies. 

Reverberation mapping (Blandford & McKee 1982) is cur¬ 
rently the most direct way to obtain BH masses in AGNs, al¬ 
though the masses are uncertain by a factor that depends on the 
poorly constrained geometry of the broad-line region (BLR). 
Reverberation mapping provides an estimate of the size of the 
BLR from the lag between the variability in the photoionizing 
continuum and the broad emission lines. With a BLR size at 
hand, one can infer a virial mass for the very central region of 
the AGN, and hence for the BH (Ho 1999; Wandel et al. 1999; 


Kaspi et al. 2000; Peterson et al. 2004): Mbh = /^blrt'/G, 
where v is the velocity dispersion of the BLR gas and the fac¬ 
tor / accounts for the geometry of the BLR (e.g., Onken et al. 
2004; Kaspi et al. 2005). Eor instance, a spherical BLR has 
an / value of 0.75 when v is measured using the full width at 
half-maximum (EWHM) of the broad line (e.g., Netzer 1990). 
Reverberation mapping has further been used to derive an em¬ 
pirical relation between AGN luminosity and BLR radius (the 
radius-luminosity relation): Rblr oc (Kaspi et al. 2000, 

2005; Greene & Ho 2005b). Since reverberation mapping data 
are time-consuming to obtain, and are currently available only 
for a small number of objects, practical applications of the virial 
technique to estimate BH masses have relied on the radius- 
luminosity relation and BLR line widths measured from single¬ 
epoch spectra (e.g., McLure & Dunlop 2001; Vestergaard 2002; 
Greene & Ho 2004). Given the many uncertainties inherent in 
the virial technique (see Krolik 2001 for a detailed discussion), 
in particular the unknown geometry of the BLR, one may justi¬ 
fiably question its reliability. On the other hand, Gebhardt et al. 
(2000b) and Eerrarese et al. (2001) showed that virial BH mass 
estimates agree surprisingly well with masses inferred from the 
Mbh-O'* relation, at least for a handful of objects. Subsequent 
work (Nelson et al. 2004; Onken et al. 2004) has increased the 
sample size, still considering a total of only 17 objects between 
the two samples. Onken et al. find that reverberation mapped 
masses, using the standard assumption of a spherical BLR and 
the EWHM, are ~ 0.26 dex below the masses inferred from 


1 


2 


GREENE & HO 


the Tremaine et al. (2002) fit to the Mbh - cr* relation; Nel¬ 
son et al. obtain essentially the same result, concluding that 
the AGN sample lie systematically below the sample of inac¬ 
tive galaxies by 0.21 dex. However, as noted by Nelson et al., 
the significance of the measured offset is low, given the large 
scatter (0.46 dex) in the data. Since this value represents the 
zeropoint in the BH mass scale for AGNs, many more cr* mea¬ 
surements are required to improve the calibration. An enlarged 
sample would allow us to resolve outstanding technical ques¬ 
tions, such as whether the EWHM or the actual second moment 
of the broad-line profile is a better measure of BLR velocity 
dispersion (e.g., Vestergaard 2002; Peterson et al. 2004), as 
yet unaddressed for virial masses obtained with single-epoch 
spectra. Eurthermore, we would be able to explore the impor¬ 
tance of BLR geometry in a statistical way, both by addressing 
the possible importance of inclination (e.g., McLure & Dun¬ 
lop 2001) and by obtaining a more secure empirical measure of 
/. The BLR geometry may even depend on properties of the 
system such as BH mass and accretion rate. With such calibra¬ 
tions in hand, we may begin to address other fundamental ques¬ 
tions, such as whether the Mbh-ct* relation varies as a function 
of BH mass (e.g., Robertson et al. 2005) or evolves with time 
(e.g.. Shields et al. 2003). 

The primary hurdle in this endeavor lies with the difficulty of 
obtaining reliable cr* measurements in active galaxies. AGNs 
are bright, and thus detectable to cosmological distances, but 
the strong continuum of Type 1 sources dilutes the starlight 
while their rich emission-line spectrum confuses and distorts 
the shape of the stellar absorption features. This paper presents 
a comprehensive discussion of how to tackle this problem. Us¬ 
ing a relatively large sample of AGNs for which we can mea¬ 
sure cr* reliably, we discuss the relative merits and complica¬ 
tions of measuring ct* from different spectral regions. The sam¬ 
ple is described in §2. Section 3 introduces the direct-fitting 
code we have developed for measuring ct* and discusses the 
battery of tests to which we have subjected it to evaluate its 
robustness and limitations. The many challenges inherent in 
dealing with AGN spectra are outlined in §4, where we devote 
considerable attention to finding the optimal spectral regions for 
measuring cr* under realistic conditions encountered in AGNs. 
We end with some practical suggestions to serve as a guide for 
other researchers (§5), followed by a summary (§6). A compan¬ 
ion paper (Greene & Ho 2005c) uses the final measurements 
obtained from this analysis to provide a new appraisal of the 
Mbh-ct* relation of AGNs. 

Throughout we assume the following cosmological param¬ 
eters to calculate distances: Hq = 100 h = 71 km s“' Mpc“\ 
Hn, = 0.27, and Ha = 0.75 (Spergel et al. 2003). 

2. THE SAMPLE 

We Utilize the large and homogeneous database provided by 
the Sloan Digital Sky Survey (SDSS; York et al. 2000), which 
will eventually obtain imaging and follow-up spectroscopy for 
about one-quarter of the sky. Briefly, spectroscopic candidates 
are selected based on multi-band imaging (Strauss et al. 2002) 
with a drift-scan camera (Gunn et al. 1998). A pair of spec¬ 
trographs is fed by 3"-diameter fibers, covering ^ 3800-9200 
A with an instrumental resolution of A/A A w 1800 (Gaussian 
iTinst ~ 70 km s“'). The spectroscopic pipeline performs ba¬ 
sic image calibrations, as well as spectral extraction, sky sub¬ 
traction, removal of atmospheric absorption bands, and wave¬ 
length and spectrophotometric calibration (Stoughton et al. 


Table 1. Sample Properties 


SDSS Name 
(1) 

z 

(2) 

S/N 

(3) 

O' 

(4) 

J000805.62-1-145023.4 

0.0454 

32 

140 

± 27 

J004236.86-104921.8 

0.0419 

44 

78.4 

± 10 

J011703.58-1-000027.3 

0.0456 

33 

98.8 

± 16 

J020615.99-001729.1 

0.0426 

46 

216 

± 30 

J033013.26-053236.0 

0.0131 

43 

99.2 

± 11 

J080243.39-1-310403.3 

0.0409 

32 

151 

± 17 

J080538.66-1-261005.4 

0.0170 

29 

100 

± 14 

J082510.23-1-375919.7 

0.0214 

30 

98.7 

± 12 

J083202.16-1-461425.7 

0.0459 

33 

104 

± 13 

J083949.64-1-484701.4 

0.0394 

38 

133 

± 12 

J092438.88-1-560746.9 

0.0250 

49 

135 

± 10 

J093812.26-1-074340.0 

0.0219 

38 

119 

O 

00 

-H 

J094529.37-1-093610.4 

0.0133 

33 

44.2 

± 15^ 

J094838.42-1-403043.7 

0.0470 

30 

91.4 

± 16 

J095033.15-1-441851.7 

0.0153 

43 

164 

± 13 

J101912.57-1-635802.7 

0.0414 

53 

239 

± 13 

J104451.73-1-063548.6 

0.0276 

33 

107 

± 9.0 

J111857.68-I-580323.5 

0.0279 

34 

50.3 

± 18^ 

J112536.16-1-542257.1 

0.0206 

37 

124 

± 12 

J120114.35-034041.0 

0.0195 

38 

50.1 

± 16^ 

J120257.81-I-045045.0 

0.0207 

18 

102 

± 20 

J121607.09-I-504930.0 

0.0308 

39 

138 

± 13 

J121754.97-1-583935.6 

0.0233 

34 

85.6 

± 14 

J122324.13-1-024044.4 

0.0236 

42 

78.1 

± 15 

J123237.48-I-662452.3 

0.0468 

48 

124 

± 9.0 

J130620.97-1-531823.1 

0.0237 

32 

90.7 

± 11 

J131305.80-1-012755.9 

0.0294 

38 

95.4 

± 13 

J134952.84-1-020445.1 

0.0328 

38 

138 

± 13 

J140018.42-1-050242.2 

0.0342 

39 

145 

± 14 

J143452.45-1-483942.7 

0.0365 

36 

138 

± 15 

J144629.97-1-500130.5 

0.0427 

25 

86.8 

± 19 

J145706.80-1-494008.4 

0.0134 

41 

108 

± 11 

J145901.35-1-611353.5 

0.0388 

34 

113 

± 13 

J150745.00-1-512710.2 

0.0458 

31 

139 

± 17 

J161156.31-1-521116.8 

0.0410 

29 

116 

± 10 

J161951.31-1-405847.2 

0.0379 

25 

89.1 

± 20 

J162012.75-1-400906.1 

0.0285 

29 



J162345.88-1-410456.3 

0.0335 

32 

35.7 

± 13^ 

J164909.58-1-361325.8 

0.0308 

36 

202 

± 22 

J232721.96-1-152437.3 

0.0457 

57 

193 

± 18 


^These measurements are below the resolution limit of 
SDSS (o'* = 70 km s~^) and thus have very low accuracy. 


Note. — Col. (1): Official SDSS name. Col. (2): Red- 
shift measured by the SDSS pipeline. Col. (3): Mean 
signal-to-noise (S/N) ratio per pixel in the Ca triplet re¬ 
gion. Col. (4): <T* (km s~^) measured from Ca triplet. 

2002). Einally, spectral classification is performed using cross¬ 
calibration with stellar, emission-line, and active galaxy tem¬ 
plates. The Third Data Release (DR3; Abazajian et al. 2005) of 
SDSS contains 45,260 spectroscopically identified AGNs with 
z<2.3. 

There are a limited number of optical stellar spectral features 
that have high enough equivalent width (EW) to be visible even 
in the presence of significant AGN contamination (see Eig. 1). 
The most common spectral regions employed are either around 
the Ca II AA8498, 8542, 8662 triplet (hereafter CaT) or the 
Mg lb AA5167, 5173, 5184 triplet. Since we plan to compare 
various spectral regions, we select all spectroscopically identi- 
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Fig. 1. — Sample spectral fits for the regions around {left column) the Ca ii AA 8498, 8542, 8662 triplet, {middle column) the Mg ib AA5167, 5173, 5184 triplet, 
and {right column) Ca H+K AA3969, 3934. The top four rows are objects from our sample; the official SDSS name is given in the upper left-hand comer for each 
row. All spectra are shown at rest wavelengths and have been normalized to their local continua. The shaded areas denote regions that were excluded from the fit. 
For each object, the top spectrum is the galaxy spectmm, the middle spectmm is our best-fit template model (offset for clarity), and the bottom spectrum is the 
residuals. Note that the CaT features are much less diluted than Mg \b. In the Mg \b region we have only fit the range 5250-5820 A, but we include the Mg \b 
features to specifically highlight their high level of AGN dilution and Fe ii contamination. In the Ca H-i-K region the entire H line needs to be masked from the fit. 
The bottom row shows the templates used in the analysis: the average stellar template, the A star template, and the Fe template derived from I Zw 1 (see text for 
details). The current Fe template does not extend to the CaT region, but Fe ii emission is known to be weak in this region for AGNs (e.g., Persson 1988). 



fied AGNs from DR3 with z < 0.05, such that at least two of the 
three CaT lines are in the band. Of the 66 objects satisfying this 
criterion, we limit our attention to those with Ca II A8542 EW 
> 1.5 A and a local signal-to-noise ratio (S/N) of 18 or higher 
per pixel. On inspection, we rejected nine candidates meeting 


these requirements because of confusion from either Paschen or 
O I A8446 emission features or strong night-sky line residuals, 
leaving a final sample of 40 objects (Table 1). 

3. STELLAR VELOCITY DISPERSION MEASUREMENTS 
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3.1. The Method 

Depending on the application, a variety of methods are used 
to measure stellar velocity dispersions. Those involving Eourier 
techniques include cross-correlation (Tonry & Davis 1979), the 
Eourier quotient (Simkin 1974; Sargent et al. 1977), and the 
Eourier correlation quotient (Bender 1990). The use of Eourier 
space is natural, as velocity broadening is reduced to a mul¬ 
tiplication, making the methods computationally cheap. It is 
also possible, however, to compare a grid of broadened tem¬ 
plate spectra directly with the galaxy spectrum in pixel space 
(Burbidge et al. 1961). While direct-fitting is comparatively 
computationally expensive, it is not particularly demanding by 
modern standards. There are many advantages to direct-fitting, 
as described, for example, by Rix & White (1992). Since the 
quality of the fit is measured in pixel space, error analysis and 
the effects of noise are easily incorporated, as is the masking 
of corrupted spectral regions or emission lines. Eurthermore, 
AGN features are easily incorporated as additional model com¬ 
ponents in a direct fit. Eollowing Barth et al. (2002a), we build 
a model, M(x), that is the convolution of a stellar template spec¬ 
trum, T(x), and a line-of-sight velocity broadening function ap¬ 
proximated as a Gaussian, G(x): 

M(x) = P{x) {[T(x) ® G(x)] + C(x)}. (1) 

The spectra and templates are normalized to the local contin¬ 
uum using identical spectral regions free of strong emission 
lines. We assume a Gaussian velocity broadening function for 
simplicity, but technically it is possible to solve for higher- 
order moments of the velocity profile (e.g., Rix & White 1992; 
van der Marel 1994). In our case the data quality does not war¬ 
rant such a treatment. C(x) is a model for the AGN continuum, 
here assumed to be a single power law, where both the nor¬ 
malization and slope are allowed to vary in the fit. This com¬ 
ponent could comprise other additive components, such as the 
Ee II emission from the BLR of the AGN that forms a “pseudo¬ 
continuum” throughout the optical spectrum (e.g., Erancis et 
al. 1991). The polynomial factor, P(x), is required to account 
for variations in continuum shape between the template and the 
galaxy (see, e.g.. Kelson et al. 2000), which, in our case, can re¬ 
sult from a combination of internal reddening in the host galaxy 
(Galactic extinction is already accounted for), differing stellar 
populations, and residual calibration errors. We typically use 
a third-order Legendre polynomial. Altogether, in addition to 
the velocity dispersion, we solve for the velocity shift, an am¬ 
plitude and slope for the power-law component, and four poly¬ 
nomial components, for a total of eight free parameters. The 
best-fit parameters are determined through a minimization of 
the statistic using the non-linear Levenberg-Marquardt min¬ 
imization algorithm as implemented by mpfit in IDL. 

Eor stellar templates, we use 32 G and K giant stars in the 
old open cluster M67 observed with SDSS (see Bernardi et al. 
2003). While the use of field stars is far more common, these 
templates are preferable because they are observed in an iden¬ 
tical manner as the program objects. The near-solar metallic- 
ity of M67 ([Ee/H] = -0.05±0.03; Montgomery et al. 1993) 
means that template mismatch problems resulting from metal- 
licity (discussed below) are comparable to those expected from 
field stars. Note that we are using a restricted range of stellar 
spectral types for our templates, while a wider range of stellar 
types must be present at some level in the galaxies. If we had 
sufficient S/N, and AGN contamination were not serious, we 
could determine the optimal combination of template stars to 



Fig. 2.— Comparison between a* measurements from SDSS (Heckman et 
al. 2004) and our own measurements of a select sample of 504 galaxies from 
that work. The agreement is reasonable, since the discrepancies at the low and 
high end correspond to regions with lower S/N spectra. The solid line repre¬ 
sents equality of the two measurements, and the dotted lines mark the nominal 
resolution limit of the SDSS. The top panel shows the fractional difference 
Sa^= [cr(ours)-f7(SDSS)]/cr(SDSS). 

better reproduce the true stellar population mix (e.g., Rix et al. 
1995; Kobulnicky & Gebhardt 2000), but such a treatment is 
beyond the scope of this work. Instead, we focus on well- 
defined spectral regions that are dominated by old stellar popu¬ 
lations in an effort to minimize template-mismatch effects. Em¬ 
pirically, as found by, for instance, Barth et al. (2003, 2005), 
early K giants consistently provide the closest match to both the 
Mg lb and CaT regions of many AGN samples. Therefore, we 
believe we are justified in this choice of stellar templates (we 
will examine potential biases in more detail in §3.2.2). One un¬ 
fortunate characteristic of these templates is their rather moder¬ 
ate S/N. While higher than the AGN spectra, they have a mean 
S/N of only 30, 75, and 65 in the regions surrounding Ca Hh-K, 
Mg lb, and CaT, respectively, resulting in additional uncertain¬ 
ties from the velocity standards themselves. We have therefore 
adopted a rather unorthodox procedure. Rather than performing 
an individual fit with each template star and then representing 
our best fit as their average (e.g., Barth et al. 2003), we cre¬ 
ate a master template by averaging all 32 stars (hereafter the 
“average template”), and use that as our velocity standard. To 
make the average template, we first align the spectra in velocity 
(shifts of ^ 0.1 pixel or 7 km s“' are typical) and then compute 
a direct average. We find from simulations that our error bars 
are decreased using this average spectrum, while no systematic 
uncertainties are introduced. Likewise, our effective resolution 
is not affected (see §3.2.2). 

3.2. Tests of the Direct-fitting Code 
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As this is the first presentation of our velocity dispersion 
code, we devote extra care to examine the reliability and limi¬ 
tations of the method. To begin with, using comparisons with 
cr* from the literature, we demonstrate that we can reproduce 
established measurements. We then perform a suite of simula¬ 
tions, in which we measure the velocity dispersion of an artifi¬ 
cially velocity-broadened input model spectrum with noise and 
A star contamination added, in order to quantify the effect of 
these variables on the accuracy and precision of our measure¬ 
ments. 

3.2.1. External Comparisons 

As an initial verification that our code functions reliably, we 
demonstrate that we can reproduce published stellar velocity 
dispersions. Heckman et al. (2004) present velocity dispersion 
measurements for SDSS Type 2 AGNs, and the measurements 
for the 33,589 AGNs from the Second Data Release are made 
publicly available by Brinchmann et al. (2004). The velocity 
dispersions were measured from the SDSS data using a direct- 
fitting algorithm to be described in detail by D. J. Schlegel et 
al. (in preparation). In short, their template stars consist of the 
first few eigenspectra from a principal component analysis of 
the echelle stellar spectra in the Elodie database (Moultaka et 
al. 2004). A Gaussian broadening function is assumed, and the 
fitting is performed over the region from 4100 to 6800 A. These 
form an excellent comparison sample for our work because they 
are based on SDSS data using a very comparable method. Note 
that because these are Type 2 AGNs, unlike our Type 1 sources, 
their continua are dominated by starlight, although there may 
be some nonstellar continuum as well (e.g., Zakamska et al. 
2003). As long as we mask the emission features, the stellar 
continuum can be fitted in a straightforward manner. In or¬ 
der to select reasonably sized subsamples of only the closest 
galaxies at each value of cr*, we applied different redshift con¬ 
straints in different cr* bins. For our chosen bins in velocity 
dispersion we employed the following (rather ad hoc) redshift 
cuts: 70-125 km s“\ z < 0.025; 125-150 km s“\ z < 0.03; 
150-170 km s-i, z < 0.035; 170-225 km s”', z < 0.04; 225- 
280 km s“', z < 0.05; and 280^00 km s“', z < 0.1. After 
removing all galaxies with a S/N < 20, we are left with a to¬ 
tal of 504 objects. We fit a slightly different spectral region, 
from 4000 to 6200 A, but we find excellent agreement with 
the SDSS measurements, as shown in Figure 2. The mean 
in the residuals, when expressed as a fractional difference, is 
((5cr*)=([cr(ours)-cr(SDSS)]/cr(SDSS)) =-0.016 ±0.13. Al¬ 
though the absolute scatter increases at the highest velocity dis¬ 
persions, in relative terms the scatter is dominated by the lower 
velocities, particularly close to the resolution limit. Dividing 
the sample, we find (5cr*) = -0.018 ±0.14 for ct* < 200 km s“^ 
versus (5cr*) = -0.013 ±0.087 at larger cr*. 

We also use our comparison with the SDSS sample to explore 
potential uncertainties introduced by the polynomial order used 
in the fit. Our standard fit employs a third-order polynomial. If 
we lower the order to 2 we find a systematic offset of roughly 
-5% from the SDSS values. Even worse, when we use a fourth- 
order polynomial we find large systematic differences. At low 
cr* (< 200 km s“') we find a ~ 25% positive offset, while there 
is reasonable agreement for cr* > 350 km s“', very similar to the 
trend reported by Barth et al. (2002a). Naively, we had expected 
the large spectral range to minimize the importance of the poly¬ 
nomial order. However, as noted by Barth et al., since our spec¬ 
tra are flux-calibrated, high-order polynomial terms may begin 


to fit 



Fig. 3.— Input-output test using an average of 10 template stars as a model 
galaxy. <5(7*= [(T(out) — (T(in)]/(7(in). The triangles are derived from the aver¬ 
age of the fit from each individual template star, while the circles come from 
fitting with the average template. The results are similar, but the average tem¬ 
plate behaves better at the lowest velocity dispersions. We see that the effect 
of resolution becomes noticeable for cr* < 100 km s“'. 

individual spectral features. To minimize this systematic effect, 
in the following we restrict our fits to polynomials of order 3. 

As an additional check, because our code is newly devel¬ 
oped for this work, we have compared our results with the 
measurements of two well-established direct-fitting codes. M. 
Sarzi kindly provided an IDL adaptation of the Rix & White 
(1992) code for intercomparison. For a subsample of 13 high- 
S/N, pure absorption-line (early-type) galaxies taken from the 
SDSS sample of Bernard! et al. (2003) we found, for 5cr* = 
[(j(ours)-(j(test)]/(T(test), that (Sa^) = 0.01 ±0.03. Since the 
Rix & White code was not designed specifically to deal with 
AGN contamination, we also compared our results with those 
obtained using the code presented in Barth et al. (2002a). A. J. 
Barth kindly used his code to fit six of the objects from our sam¬ 
ple over an identical fitting region. Again, our results agreed 
within the uncertainties, with no systematic differences appar¬ 
ent, yielding {6a^,) = 0.01 ±0.05, where we attribute the small 
differences to differences in the minimization routines. These 
tests give us additional confidence in the reliability of our code. 

3.2.2. Systematics 

Many factors may contribute to uncertainties in our ct* mea¬ 
surements, including template mismatch, resolution, and S/N 
(we discuss emission from the AGN separately). By generat¬ 
ing model spectra, we can examine how each of these factors in 
turn impacts our ability to recover the input velocity width. We 
build model galaxies as linear combinations of template spec¬ 
tra, which we then broaden to a range of widths between a = 
50 and 400 km s“', assuming a Gaussian line-of-sight veloc¬ 
ity distribution. We then vary parameters of interest, such as 
S/N, and evaluate our ability to recover cr*, using as the figure 
of merit the fractional error (5 ct* = [(T(out)-cr(in)]/cr(in). For 
later use, we focus on (1) the region around Ca Hh-K (3910- 
4060 A) with Ca H (3955-3985 A) masked, (2) the region red- 
ward of Mg lb (5250-5820 A; hereafter called the “Fe region”) 





/ km s"*- 



100 200 300 400 

/ km 


Fig. 4.— Uncertainty, expressed as a fractional difference 5cr*= [cr(out) — cr(in)]/cr(in), due to template mismatch for (a) CaT, {b) the Fe region (5250-5820 A), 
and (c) Ca K. We create 100 different model galaxies from random collections of 10 template stars, and depict one standard deviation above and below the mean 
5a* by the shaded region. The mean is demarcated by the thin solid line. We also run simulations with 10% {circles), 20% {squares), and 30% {triangles) A star 
contribution. The points represent the mean 5a*, and the error bars are one standard deviation in 5a*. Because of the strange behavior of Fe shown in {b), we 
present alternate simulations in {d) excluding the six largest outliers. 


with [Fe VII] A5721 masked, and (3) the region around the CaT 
(8470-8700 A). Each fit is done locally, without attempting to 
match the polynomial or AGN power-law between regions. 

To begin with the simplest possible test, we fit a model 
galaxy, composed of 10 template stars with only artificial 
broadening, in the CaT region. We are using a single average 
template in order to increase the effective S/N of our template, 
but for comparison we also compute ct* as the average of the 
values obtained from each of the 32 individual template stars. 
The resulting (icr* values for both approaches are shown in Fig¬ 
ure 3. No apparent offset is introduced by the use of an average 
template, while the fractional error close to the resolution limit 
is diminished. In general, while additional uncertainty is intro¬ 
duced below cr* w 100 km s“' due to resolution, it is only on 
the order of ~ 2%. We follow previous work (e.g., Bernardi et 


al. 2003; Heckman et al. 2004) and limit our attention only to 
velocity dispersions larger than the nominal instrumental reso¬ 
lution of SDSS, which we take to be crinst = 70 km s“'. 

Even for predominantly old stellar populations (G and K gi¬ 
ant stars), template mismatch will introduce feature-dependent 
uncertainties. To investigate the magnitude of these uncertain¬ 
ties for each fitting region, we generate a family of 100 model 
galaxies. As above, the model consists of 10 templates, which 
are selected randomly for each run, and then coadded and artifi¬ 
cially broadened. The shaded region in Figure 4 represents one 
standard deviation around Jct* for the ensemble of 100 trials. 
Note that the CaT region is well-behaved (Fig. Aa), while the 
Fe region is, on average, ~ 2% low (Fig. Ab). However, when 
we remove the six largest outlying templates from the parent set 
of 32 and repeat the simulations, the negative offset is removed 
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Fig. 5.— Input-output simulations of velocity dispersion measurements from (a) the CaT region, (b) Fe (5250-5820 A), and (c) Ca K of galaxy models con¬ 
structed as the average of 10 template stars, = [cr(out)-cr(in)]/cr(in). Simulations are shown for S/N of 15 to 55, in steps of 5. The highest S/N (i.e., no 
additional noise) is ~ 205, 240, and 95 for CaT, Fe, and Ca K, respectively (shown in open circles). We highlight in solid lines the cases for S/N = 30 {triangles) 
and S/N = 50 {squares), since the majority of our data are spanned by these limits. The middle and right columns correspond to /agn = 50% and 66%, respectively. 


(Fig. Ad). In the Ca H+K region Jct* has a mean value of zero, 
but by far the largest dispersion. 

While this experiment gives us a feel for template mismatch 
among old stellar populations, additional uncertainty is intro¬ 
duced by young stellar populations. Particularly, it has been 
shown that the host galaxies of luminous AGNs do not con¬ 
tain a uniformly old population, but rather are more likely to 
contain a post-starburst population than quiescent galaxies of 
similar total stellar mass (e.g., Terlevich et al. 1990; Kauff- 
mann et al. 2003b; Nelson et al. 2004). Due to rotational and 


pressure broadening, a significant A star population can sig¬ 
nificantly broaden the hydrogen lines, and potentially bias the 
measurements in the regions around CaT and Ca Hh-K. There¬ 
fore, we run an additional set of simulations in which we ex¬ 
plicitly include an A star contribution, using an A star spectrum 
observed by SDSS, Gaussian-broadened to the same width. In 
Figure 1 (bottom panel) we show that the A star spectrum has 
significant features in both the Ca Hh-K and CaT regions, while 
it is nearly featureless in the Fe region. In the case of CaT, 
the features are a blend of Paschen and CaT lines, while the 
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Ca K line is apparent to the left of a blend of Ca H and He in 
the Ca H+K region. The (5 ct* values for models including A 
star contributions of 10%, 20%, and 30% are shown in Eigure 
4. As expected, both the Ca H+K and CaT regions are more 
sensitive to the presence of A stars, although the impact is rel¬ 
atively small (< 10%). Note, however, that the induced spread 
is largest in the case of Ca H+K. 

The situation is more puzzling in the case of the Ee region, 
where the uncertainties appear to decrease with the addition of 
an A star component. Apparently an A star contribution leads 
to a better fit of the outlying templates, which were removed 
to generate Eigure Ad. Various scenarios may account for this 
behavior. If, for instance, the outliers have a continuum shape 
opposite to that of an A star (i.e., very red), then the addition 
of an A star would cancel them out. We investigated this possi¬ 
bility by running simulations including the A star but excluding 
the outlying templates, shown as points in Eigure Ad, but the in¬ 
duced errors do not increase, so this is probably not the solution. 
Another possibility is that the outlying templates have stellar 
absorption features with larger EWs, so that the additional dilu¬ 
tion from the A star continuum improves their fit with the rest of 
the templates. We measured the EW in Mg \b (using the Lick 
index definition; e.g., Trager et al. 1998) for each individual 
template, and the mean EW for the outliers is 5.2 A, compared 
to (EW) = 3.3 A for the rest, indicating this is the most likely 
explanation. Since our actual galaxies will probably contain 
some A star contribution, we have decided to retain all 32 tem¬ 
plates in our average template, despite the apparent skewness 
of Eigure Ab. We verified that our actual fits are not strongly af¬ 
fected if we use an average template excluding the six outliers. 
Between the nominal results and those derived from an average 
template excluding the six outliers, we find an average relative 
difference of ((5cr*) =-0.06 ±0.1. 

We now turn to the impact of varying S/N on the reliability 
of our measurements. Eor these simulations, we begin with a 
well-behaved model galaxy from the simulation above. We de¬ 
grade the S/N by adding mean-zero Gaussian random deviates, 
whose uncertainty at each pixel is given by the SDSS pipeline 
error array, and whose amplitude is varied to decrease the ef¬ 
fective S/N from 60 to 15, in steps of 5. Note that the model 
galaxy has a S/N of ~ 205, 240, and 95 for the CaT, Ee, and 
Ca H+K regions, respectively, which forms the upper limit on 
our noise simulations. In practice, the majority of our objects 
are found between a S/N of 30 and 50 (although the distribu¬ 
tion is somewhat lower in the Ca H+K region). Eor each S/N, 
we generate 50 model spectra, in order to average over the vari¬ 
ations caused, for instance, when a noise spike falls on a crucial 
absorption feature. We do not find a systematic bias in Jct* from 
variations in S/N or dilution and so for simplicity we report the 
standard deviation in 5(T*over the 50 trials, denoted RMS(()ct*). 
The results are shown in the left-hand panels of Eigure 5, for 
each spectral region of interest. Erom these simulations it is 
clear that even at the modest S/N of 20, which is common for 
SDSS spectra, it is possible to achieve ct* measurements with 
~ 20% accuracy. 


4. MEASURING cr* IN AGNS 

Inspired recently by the discovery of the Mbh - ct* relation, 
many groups have measured cr* in Type 1 AGNs (e.g., Jimenez- 
Benito et al. 2000; Eerrarese et al. 2001; Barth et al. 2002b, 


2003, 2005; Nelson et al. 2004; Onken et al. 2004; Treu et 
al. 2004; Woo et al. 2004; Botte et al. 2005; Garcia-Rissmann 
et al. 2005). As noted by many of these authors, measuring 
cr* in the presence of an AGN is particularly challenging. The 
AGN power-law continuum dilutes the stellar features, effec¬ 
tively lowering the S/N, while emission lines from the AGN can 
add (sometimes subtle) biases into the fits. In order to mitigate 
these problems, the above studies fit specific spectral regions 
with minimal AGN contamination but high-EW stellar absorp¬ 
tion features, rather than attempting to find global solutions for 
large regions of the spectrum. We adopt the same approach. In 
particular, the utility of the CaT region has long been recog¬ 
nized (Pritchet 1978; Dressier 1984). The CaT lines are intrin¬ 
sically narrow compared to galactic velocity dispersions, with 
intrinsic widths of cr w 20-25 km s“' for K giant stars, ranging 
up to ~ 30 km s“' for supergiant stars (Eilippenko & Ho 2003; 
Martini & Ho 2004). Eurthermore, the AGN continuum is min¬ 
imized in their vicinity, since AGN continua are typically blue. 
Although higher-order Paschen lines (P13 A8503, P15 A8545, 
and PI6 A8665) overlap with the CaT lines, these features are 
rarely seen in emission (an exception is NGC 4395, which, as 
noted by Eilippenko & Ho, has anomalously high-EW emission 
features). Emission from CaT itself can fill in the absorption 
features as well, but these objects were eliminated by selection. 
Einally, there may also be OI A8446 or [Ee II] A8616 emission 
from the AGN. In principle we can mask either of these features 
in a straightforward manner within our fit, but again, these lines 
are weak by selection in our sample. We measured meaning¬ 
ful cr* measurements for all but one of our objects, although an 
additional four are below our nominal resolution limit of ^ 70 
km s“^ (Table 1). At the S/N of the SDSS data, our ability to 
recover accurate cr* measurements below the resolution limit is 
significantly decreased, and so we conservatively do not include 
these four measurements in our analysis. 

Of all high-EW features in the optical bandpass, the CaT 
lines are simultaneously the least affected by template mis¬ 
match and AGN contamination, making them our “gold stan¬ 
dard.” Unfortunately, because they occur in the red region 
of the spectrum, they are inaccessible to optical bandpasses 
above z « 0.05. Eor this reason, we must examine the utility 
of other spectral regions with reasonably strong absorption fea¬ 
tures. We will focus on two: the region surrounding the Mg \b 
triplet (5040-5820 A) and the region around the Ca H+K lines 
(3900^060 A). We measure cr* using these alternate spectral 
regions, and compare the results directly with the CaT val¬ 
ues. In order to evaluate whether a given test region matches 
CaT, we use the same figure of merit as defined above, (5cr* = 
[cr(test)-cr(CaT)]/cr(CaT), using in all cases the 35 measure¬ 
ments with cr* (CaT) > 70 km s“'. 

4.1. Power-law Continuum Contamination 

Relatively speaking, it is straightforward to model the con¬ 
taminating effect of a pure power-law component on the final 
cr* measurement. We perform simulations identical to those 
described in §3.2.2, starting with a model composed of 10 tem¬ 
plate stars, adjusting the S/N from 15 to 60, in steps of 5, and 
measuring Jcr*, except that now we also include a power-law 
component to represent the AGN continuum. We use a single 
power law, oc v~^, comprising 0, 33, 50, 60, 66, 70, 75, 80% 
of the total local continuum, as measured in the 20 A at each 
end of the band. Assuming our nominal power-law slope of-1, 
an AGN fraction of 50% (66%) in the Ee region corresponds to 
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Fig. 6.— Examples of the Mg \b spectral region for ob¬ 
jects with significant Fe ii contamination, as well as, in some cases, 
contamination from [Fe vii] A5158, [Fe vi] A5176, [N i] AA5197, 
5177, [Fe vii] A5278, [Fe xiv] A5303, and [Fe vi] A5335. Plot¬ 
ted in order of decreasing [Fe vii] strength, from top to bottom, are 
SDSS J162012.75-I-400906.1, J120114.35-034041.0, J083949.65-H484701.4, 
J112536.16-I-542257.1, J150745.00-I-512710.2, and J101912.56-H635802.6. 

48% (65%) AGN fraction in the CaT region and 73% (84%) 
AGN fraction in the Ca Hh-K region. For general comparison 
with the literature, we also compute the EW in Mg \b, CaT, and 
Ca Hh-K from the average template star for AGN fractions of 
0%,50%, and 66% in the Fe region. We use the Lick index def¬ 
inition of the Mg \b index (Trager et al. 1998), the Terlevich et 
al. (1989) definition of the CaT A8542 A line, and the Brodie 
& Hanes (1996) definition of the HK index. For no dilution we 
find EW = 4, 3, and 22 A for Mg lb, CaT, and Ca Hh-K, respec¬ 
tively. These decrease to EW = 2, 1, and 6 A for /agn = 50%, 
and to EW = 1, 0.9, and 4 A for /agn = 66%. We show repre¬ 
sentative cases with 50% and 66% contamination in the middle 
and right panels of Figure 5. Although we have run a larger 
suite of simulations, these examples bracket the majority of the 
measured values for our sample. We assign an uncertainty to 
each of our galaxies by interpolating between the grid points 
from the simulations to its measured S/N and dilution. These 
simulations show that if template mismatch and AGN emission 
features were not present, we would achieve the best measure¬ 
ments using the line with the highest EW (namely Ca Hh-K in 
this case). While not directly relevant to this particular sample, 
it is of general interest to determine the limiting AGN dilution 
for each spectral feature. As a general indicator, we set the max¬ 
imum AGN dilution at the point where 5cr* = 1 at a S/N = 20. 
This occurs at an AGN fraction of 71 % for CaT, 85% for the Fe 
region, and 90% for Ca Hh-K. Although the Ca Hh-K lines en¬ 
dure the highest AGN contamination, they are very susceptible 
to systematic uncertainties induced by template mismatch, and 
so we prefer the CaT lines despite their lower EWs. 

As a sanity check, we compare to the BE Lac objects studied 
by Barth et al. (2003), whose dominant sources of uncertainty 
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are S/N and dilution from a pure power-law continuum. For 
their sample, with dilutions of 50%-80% and S/N between 50 
and 300, they find ~2%-ll% errors using either the CaT or 
Mg lb spectral regions. Reassuringly, to the extent that they 
overlap, our simulations are consistent with their results. 

4.2. Emission-line Contamination 

As we have seen, contamination from a power-law contin¬ 
uum only effectively lowers the S/N of the measurement. Far 
more pernicious are the subtle biases that may be introduced 
by AGN emission features and template mismatch. Between 
the Mg lb lines themselves and the prominent Fe I blends at 
5270 A and 5335 A, this region contains sufficiently high-EW 
features for measuring velocity dispersion, but they are sus¬ 
ceptible to template mismatch due to metallicity effects. Be¬ 
cause of the well-known correlation between [Mg/Fe] and ct* 
(e.g., O’Connell 1976; Worthey et al. 1992; Trager et al. 1998; 
Jprgensen 1999; Kuntschner et al. 2001), it becomes increas¬ 
ingly difficult to simultaneously fit both the Mg lb and Fe fea¬ 
tures using template spectra with solar metallicities and abun¬ 
dance ratios. The errors incurred are small, but increase with 
increasing cr*, tending to artificially increase the inferred ve¬ 
locity dispersion (Barth et al. 2002a). 

Perhaps more relevant to this study, there is significant AGN 
line contamination over this wavelength range. There are a 
number of narrow emission lines, including [Fe VII] A5158, 
[Fe VI] A5176, [N I] AA5197, 5200, [Fe VII] A5278, 
[Fe XIV] A5303, and [Fe VI] A5335 (see Fig. 6). The first 
three lines are especially troublesome, occurring as they do 
amidst the Mg lb features, and they are present in ^ 50% of 
our sample. In addition to narrow features, there is a broad 
Fe II “pseudo-continuum,” extending from ~5050 to 5520 A 
(e.g., Francis et al. 1991). When the Fe II emission is broad 
enough to be smooth, it simply lowers the effective S/N of 
the spectrum without introducing significant systematic errors 
(e.g.. Nelson & Whittle 1995). However, when the broad-line 
velocity is low enough that the individual Fe II components are 
relatively resolved, then the resulting ct* measurements may 
be biased. Some examples of galaxies with significant broad 
Fe II contamination as well as strong coronal emission lines are 
shown in Figure 6. 

Following Barth et al. (2002a), we use the region extend¬ 
ing 5040 to 5430 A, and we mask the [N I] feature from 
5190-5210 A. In Figure la, we compare the velocity disper¬ 
sions obtained from this fit with those from CaT. As is im¬ 
mediately clear, CT(Mg lb) is significantly biased. In particu¬ 
lar, there is a plume of objects with CT(CaT) values between 
75 and 150 km s“' and a(Mglb) w 80 km s“'. We find = 
([cr(MgI/7)-CT(CaT)]/CT(CaT)} = -0.23 ± 0.32. A combination 
of broad, permitted Fe II emission and high-ionization forbid¬ 
den emission must be responsible for this dramatic offset, since 
mismatches in [Mg/Fe] tend to go in the opposite direction. We 
use simulations to gauge the systematic effects due to Fe II con¬ 
tamination, and we experiment with different spectral regions to 
minimize the effects of coronal emission. 

We examine the role of broad Fe II contamination by running 
simulations, similar to those presented in §3.2.2. In the usual 
way, we create a model galaxy using 10 of our template stars 
and Gaussian-broaden it artificially to widths between 50 and 
400 km s“'. As a model of the Fe II contamination we utilize 
the Fe II spectrum of I Zw 1 (Boroson & Green 1992), kindly 
provided by T. A. Boroson. I Zw 1 is a narrow-line Seyfert 
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Fig. 7.— A comparison between cr* from the CaT region and (a) the Mg ib region (5040-5430 A) and {b) the Fe region (5250-5820 A). Note the large systematic 
offset at low cr* in the Mg ib region. The solid line represents equality, while the dotted lines mark the nominal resolution limit of the SDSS. 


1 galaxy, which for our purposes ensures that the Ee II emis¬ 
sion is both strong and narrower than the Ee II emission in our 
objects. The challenge is simulating a representative range of 
EWs and velocity dispersions for the Ee II template that accu¬ 
rately reflect what is found in real AGNs. In the sources of 
interest the Mg lb features are strong (and thus measurable), 
making it difficult to accurately determine the Ee II EW, at least 
in a model-independent way. We have chosen to vary the Ee II 
amplitude between an EW of 18 and 54 A, comparable to what 
is seen in relatively weak Ee II emitters. Eor instance, Boroson 
& Green (1992) find a range of Ee II EWs between 0 and 114 
A for their sample of nearby quasars. As demonstrated in Eig- 
ure 8, this level of Ee II contamination is hard to discern based 
on the appearance of the absorption features. However, despite 
the low amplitude, this quantity of Ee II contamination can se¬ 
riously bias the fits. Ear more important than the amplitude of 
contamination is the width of the Ee II lines. When they are nar¬ 
row enough to have significant structure, they most effectively 
confuse the fitting procedure. Empirically, the width of the Ee II 
emission is found to match that of the broad component of H/3 
(e.g., Boroson & Green 1992). Therefore, using the Mbh-ct* 
relation, for any given cr* the range of possible Ee II widths is set 
by the luminosity of the AGN through the line width-luminosity 
relation (Mbh cjc L° ®"^EWHMh/ 3 ^; Kaspi et al. 2005; Greene & 
Ho 2005b). Since, for a given Mbh the maximum luminos¬ 
ity is set by the Eddington limit, it is convenient to parameter¬ 
ize the AGN luminosity using the Eddington ratio, Lboi/TEdd, 
where Ledd = 1-26 x 10^* (Mbh/Mq) ergs s“^. In practice, the 
bolometric luminosity is difficult to obtain, and so we assume a 
canonical bolometric correction of 0.1 between the optical and 
bolometric luminosities (e.g., Elvis et al. 1994). Eor a fixed cr* 
(or equivalently Mbh), as the luminosity of the system increases 


the line width must decrease, such that the narrowest Ee II lines 
are found at the highest Lboi/Tedd- Hence, the worst contamina¬ 
tion occurs at the Eddington limit, as shown in Eigure 9, where 
we consider Lboi/TEdd = 0.01, 0.1, and 1. Also, for any given 
value of Lboi/TEdd, the bias is worst near the resolution limit.' 
We do not show the simulations with Ee II EW of 18 or 54 A, 
but the general trends are the same, with a maximum bias in the 
Mg lb region of ~ 10% for an EW of 18 A and ^ 30% for an 
EW of 54 A. In addition to the 5040-5430 A region, we also 
consider the region 5250-5820 A, which eliminates the Mg lb 
features themselves but includes the Ee I absorption features. 
As we will show below, this is a better spectral region to use in 
cases with strong [N I] and coronal emission from [Ee VII]. We 
find that this spectral region is somewhat better behaved at low 
(T* and low Eddington ratio, but worse at large cr* and at high 
Tboi/Tudd (Fig. 9b). Finally, we consider the effects of Fe II 
contamination on the Ca Hh-K lines. The Fe II features near 
~ 3950 A are both weaker and broader (see bottom of Fig. 1), 
which means that the maximum incurred bias is < 20% (Fig. 
9c; ~ 10% for EW of 18 A and - 30% for an EW of 54 A), ex¬ 
cept near the resolution limit. We note that while we have run a 
suite of simulations including many combinations of Lboi/TEdd 
and EW, in nature the Fe II strength is found to increase with 
Tboi/Tudd (e-g-, Boroson 2002), further contaminating sources 
close to their Eddington limit. 

Given the rather large bias that Fe II emission can apparently 
cause, especially at high Eddington ratios, we seek a general 
method to salvage the Mg lb region by including an Fe II model 
directly in our direct-fitting code. Equation 1, describing our 
fitting procedure, includes a term C(x), which is a model of the 
AGN continuum. Previously we have modeled the AGN com¬ 
ponent as a pure power law, but now we also include the I Zw 1 


’This explains why, for instance, Barth et al. (2005) see no bias between their Mgib and CaT measurements, since they are measuring velocity dispersions 3-5 
times the effective spectral resolution. 
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Fig. 8.— Examples of simulations including Fe ii contamination. The stel¬ 
lar continuum is modeled as the average of 10 template stars, while the Fe ii 
contamination is modeled using the Fe ii template derived from 1 Zw 1. The 
template is scaled to Fe ii EWs of 18, 36, and 54 A in our models, while the 
line width is parameterized using the Eddington ratio (Lbol/f-Edci; see discus¬ 
sion in text). From top to bottom we illustrate Lbol/f-Edd = 0.01, 0.1, and 1 with 
EW(Fe ii) = 36 A, and f-bol/^Edd = 1 with EW(Fe ii) = 54 A. The bottom spec¬ 
trum represents the most extreme contamination examined in our simulations. 
Note that in all cases the effects of the Fe ii on the stellar features are subtle 
(for instance, there is a change in the relative strengths of the Mg ib triplet). 

Fe II template. In so doing we introduce three new degrees of 
freedom: a velocity shift, a velocity dispersion, and an ampli¬ 
tude for the Fe II spectrum. Rather than allowing these param¬ 
eters to vary freely, we adopt an iterative process. Initially we 
perform the standard fit as described above with no Fe II com¬ 
ponent. When Fe II contamination is present, the residuals of 
the initial fit contain significant Fe II features. We perform a 
least-squares fit between the residuals and the I Zw 1 template 
to measure the velocity shift of Fe II, using the region from 5040 
to 5820 A to improve the leverage on the velocity. The shift is 
then constrained to this value in the final fit, while the Fe II 
width is fixed to the velocity dispersion of H/3^ (when we allow 
this parameter to be free, the code often fails). In simulations 
in which we artificially add Fe II contamination and then use 
this multi-step fitting procedure to remove it, we find that we 
are able to recover the input cr* values, with a slight negative 
offset of (5cr*w -5%, for both spectral regions. This, however, 
assumes that we know the width of the Fe II lines perfectly. 
Once we introduce a 20% variation in the Fe II width, the out¬ 
put dispersions have associated errors of 10%. 


On the other hand, when we attempt to use this fitting pro¬ 
cedure with our actual data, the results are not very encour¬ 
aging. A slight positive systematic difference of 2%-5% is 
introduced between the fits with and without the Fe II mod¬ 
eling. As we have noted, our objects are not high-Lboi/TEdd 
objects. Recall that we selected objects with both high CaT 
absorption EWs and neglible CaT emission. Since CaT emis¬ 
sion is correlated with Fe II emission, which in turn is corre¬ 
lated with Lboi/TEdd (Persson 1988; Persson & Ferland 1989; 
Boroson 2002), we have predominantly selected AGNs with 
low Fe II EW and low Lboi/TEdd- Using virial mass estimates 
based on Ha line width and luminosity (Greene & Ho 2005b), 
we find (log Lboi/TEdd) = “1-3 ±0.4 for this sample (Greene 
& Ho 2005c). In this regime, when the Ee II emission is very 
smooth, it is somewhat degenerate with the polynomial in the 
fit. While the Ee II EW estimates from the code range from 1- 
300 A, with a median value of 80 A, we do not believe these 
values are well-constrained. It is possible that in more extreme 
cases, in which the Ee II emission is more pronounced, it may 
be possible to obtain a more robust fit for the Ee II emission, and 
thus recover . Of course, high-resolution and high-S/N ob¬ 
servations would also aid significantly in performing this pro¬ 
cedure. 

We conclude that, while broad Ee II contamination can cause 
large biases in Mg lb measurements, it is not responsible for 
the offsets observed in Eigure la between Mg lb and CaT. In¬ 
stead, narrow coronal emission is probably the primary culprit, 
which suggests that alternate spectral regions with less coronal 
emission contamination and, ideally, less Ee II emission as well, 
need to be explored. The optical Ee II multiplets become neg¬ 
ligible at ~ 5520 A, and the region redward of Mg lb is rich in 
metal lines, which, while having lower EWs, in principle may 
be used to measure ct*. We have experimented with the follow¬ 
ing spectral regions: 5520-5820 A and 5520-6280 A, which 
avoid Ee II contamination almost completely, as well as both 
5040-5820 A and 5040-6280 A. Einally, we tried a compro¬ 
mise region, 5250-5820 A, which contains the high-EW Ee I 
absorption features, at the expense of some Ee II contamination. 
In all cases, we mask 5710-5740 A to remove the often-present 
[Ee VII] A5721 line, and, where relevant, we mask 5820-5950 
A to exclude the Na D feature, which is susceptible to strong 
interstellar contamination (e.g., Bica et al. 1991). In cases in¬ 
cluding the Mg lb and Ee II features we have experimented with 
masking the Mg lb features themselves, as well as other regions 
most contaminated by Ee II emission (specifically 5250-5325 
A). We have found that including the spectral region redward 
of Na D introduces significant template mismatch problems, 
and therefore do not consider it further. Of all the possibili¬ 
ties remaining, we find that the region 5250-5820 A minimizes 
the systematic offset with CaT (as shown in Eig. lb). With 
((5cr*) = ([CT(Ee)-CT(CaT)]/CT(CaT)) = -0.018± 0.20, this re¬ 
gion is unambiguously superior to Mg lb both in terms of sys¬ 
tematic offset and total scatter. Because narrow coronal emis¬ 
sion is common and introduces a strong systematic uncertainty, 
we recommend using the Ee region rather than the Mg lb region 
itself. However, we note that, at least with the moderate reso¬ 
lution and S/N of SDSS data, it is preferable to avoid both the 
Ee and Mg lb regions when either the AGN fraction is ^85% 


^We must constrain the Fe ii velocity because large velocity shifts have been observed between the BLR and the systemic velocity (e.g., Marziani et al. 1996). 
To fit the H/3 velocity dispersion, we first model and remove the narrow H/3 component using a fit constrained by the narrow Ha line fit, as described in Greene 
& Ho (2005b). The [O iii] AA4959, 5007 lines are simultaneously modeled and removed. We then fit as many Gaussian components as needed to model the broad 
component, and measure the line width from the broad-line model. 
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Fig. 9.— The effect of Fe ii contamination on cr*. We compare the systematic offsets introduced in three different spectral regions: (a) Mg ib, (b) the Fe region 
redward of Mg ib, and (c) Ca K. The model spectrum used in the simulations is the average of 10 template stars, with Fe ii contamination added using the I Zw 1 
template. The amplitude of Fe ii is kept constant for these simulations with an EW of 36 A, but its velocity width decreases as the Eddington ratio increases. We use 
the input cr* value to calculate M^u, and the assumed Eboj/Lndd value sets the continuum luminosity. We then infer FWHM of H/3, which empirically matches the 
FWHM of Fe ii (e.g., Boroson & Green 1992). The percentage systematic error introduced is quantified by Scr* = [cr(out) — cr(in)]/(j(in). 


or Lbol/r^Edd ~ 1- 

4.3. CaH+K 

The Ca H+K lines are also high-EW spectral features poten¬ 
tially useful for obtaining ct* measurements. Traditionally they 
have been avoided for a variety of reasons. Perhaps the most 
serious complication is that their line shapes are a strong func¬ 
tion of spectral type, causing systematic errors from template 
mismatch (e.g., Kormendy & Illingworth 1982; Kobulnicky & 
Gebhardt 2000). Since He overlaps with Ca H and H8 sits di¬ 
rectly blueward of Ca K, young stellar populations may seri¬ 
ously affect the region (e.g., Gebhardt et al. 2003). Apart from 
line shape, there is a very steep gradient in the local continuum 
(the “4000 A break”) whose slope depends on spectral type. An 
additional problem is caused by possible interstellar absorption 
in the host galaxy, which would lead to narrow cores (see top 
panel in Eig. 1). Einally, the lines have limited utility at low 
cr* because they are, relatively speaking, intrinsically broad. In 
combination, these factors cause the Ca Hh-K lines to system¬ 
atically overestimate ct*, as found by both Kormendy & Illing¬ 
worth (1982) and Bernard! et al. (2003). Dressier (1979), on 
the other hand, finds no systematic offsets between the Ca Hh-K 
region and the G-band (itself sensitive to template mismatch). 
However, Dressier was observing a cD galaxy, which, due to 
the (presumably exclusively) old stellar population and large 
velocity dispersion (> 375 km s“') offers the most favorable 
conditions for using Ca Hh-K. Eor a sample of late-type galax¬ 
ies, Kobulnicky & Gebhardt (2000) obtain reasonable agree¬ 
ment (~ 20%) between Ca Hh-K absorption line widths and 
kinematics inferred from H I and [O II] A3727 emission line 
widths. Note, however, that a comparison with emission fea¬ 
tures is not the same as a direct comparison with stellar absorp¬ 
tion features. Also, as emphasized by Kobulnicky & Gebhardt, 
their success hinges on their ability to include a wide range of 
spectral types in the modeling. 

^Of course, there are less direct substitutions for (j*, such as [O iii] line width 
have large scatter (see, e.g., Greene & Ho 2005a). 


Of course, AGN contamination only makes the situation 
more complicated. The power-law continuum is rising steeply 
to the blue, fractionally increasing the AGN contamination 
there, and He (and in principle [Ne III] A3968) emission ren¬ 
ders the Ca H line completely unusable. On the positive side, 
compared to the Mg lb region, as we saw in §4.2, Ee II contami¬ 
nation is relatively minimal. Despite these many complications, 
Ca Hh-K presents two major benefits compared to all other spec¬ 
tral regions. Eirst of all, it is the only stellar feature with suf¬ 
ficient EW to persist in high-luminosity AGNs; for instance, it 
is detected in the composite SDSS quasar spectrum of Vanden 
Berk et al. (2001). Secondly, the lines are blue enough to re¬ 
main in the optical bandpass beyond a redshift of 1. Ca Hh-K 
represents our only hope of obtaining cr* measurements in pow¬ 
erful AGNs or those at intermediate redshift^. Below we inves¬ 
tigate whether, despite its problems, the Ca Hh-K region can 
provide a useful diagnostic of cr* in AGNs. 

We follow the approach of the previous section, and compare 
cr* from Ca Hh-K with that from CaT to seek the combination 
of spectral region, polynomial order, and masking that yields 
the best agreement. In terms of the best spectral fitting region, 
we have found that broad He emission is always a contaminant 
of Ca H. Even when the emission line is not clearly visible, 
the centroid of the Ca H line is clearly displaced from its ex¬ 
pected position relative to Ca K, so that the code cannot find a 
reasonable solution. Eor that reason, we simply mask the Ca 
H line (3955-3985 A) in all fits; hereafter, our discussion will 
focus exclusively on the Ca K line. Over this very restricted 
spectral region, AGN dilution can be somewhat degenerate with 
increasing velocity dispersion. To mitigate this problem, we in¬ 
clude 50 A of continuum redward of Ca H in order to help the 
code determine the general level of AGN contamination. Be¬ 
cause of the limited spectral region, we also experimented with 
fewer polynomial orders or a power-law continuum constant 

(e.g., Shields et al. 2003). While statistically these techniques are accurate, they 
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Fig. 10.— (a) A comparison of cr* using the CaT and Ca K spectral regions. The solid line represents perfect agreement, while the dotted lines mark the nominal 
resolution limit of the SDSS. Error bars are derived as described in §3.2.2. (b) Same as (a), except that we have corrected the Ca K measurements using Equation 2. 
It is clear that the correction is too large at low cr*. 


with wavelength, but in both cases the final values were el¬ 
evated, and the fits were unsatisfactory. Finally, we attempted 
a two-step fitting approach to try to mitigate the potential bias 
caused by a narrow core from interstellar absorption. The first 
iteration was as described above, while for the second iteration 
we masked the line core but simultaneously fixed the velocity 
shift. Unfortunately, when only the wings are available, the cr* 
measurements become extremely (unphysically) broad, so this 
technique is untenable. In the end, we settled on a fitting re¬ 
gion of 3910^060 A with Ca H masked. The results of the 
comparison with CaT are shown in Figure 10. We find ((5 ct*) = 
([cr(Ca K)-cr(CaT)]/cr(CaT)) = 0.049 ± 0.29 km s”', which is 
fairly reasonable agreement. 

As noted above, the Ca Hh-K lines have been known to over¬ 
estimate cr*, even in the absence of AGN contamination. To 
evaluate the type of systematic errors that may ensue over a 
larger range of cr* in the absence of AGN contamination, we 
measured tT(Ca K) for the sample of 504 SDSS Type 2 Seyfert 
galaxies described in §3.2.1. As can be seen in Figure 11, 
(j(Ca K) is systematically biased, in a (T*-dependent way. At 
low cr* the Ca K values are too large, while at large values of 
(j* the Ca K values are too small, with the change occurring at 
(j* « 175kms“^. Since the comparison for our program objects 
is between CaT and Ca K, it may be more appropriate to investi¬ 
gate offsets between these two spectral regions. We have there¬ 
fore compiled an alternate sample of 374 objects from Heck¬ 
man et al. (2004) that have z < 0.03 and a S/N >10 around the 
CaT. Unfortunately, because of the redshift constraint, the ve¬ 
locity dispersions of the sample are low (< 225 km s“^), and so 
we cannot investigate global behavior from these points alone, 
but the qualitative offset is similar to that seen in the large sam¬ 
ple (Fig. 12). This trend is somewhat reminiscent of the bias 
apparent in (T*(Ca Hh-K) shown in Figure 4c due to an A star 
component, which suggests that the bias may be due to 



Fig. 11.— Comparison between tj* from 4000-6200 A as compared to 
Ca K, for a sample of 504 Type 2 Seyfert galaxies from Heckman et al. (2004). 
The solid line represents equality of the two measurements, and the dotted lines 
mark the nominal SDSS resolution limit of 70 km . The dashed line shows 
our linear least-squares fit (Equation 2). The top panel plots the fractional 
difference Sa^= [cr(Ca K)-cr(4000-6200)]/cr(4000-6200). 

differing A star contaminations as a function of a^. We inves¬ 
tigate whether young stellar populations are responsible for the 
observed offsets by seeking a correlation between Jan. and the 
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C7^ (CaT) / km (CaT) / km s-‘ 

Fig. 12.— (a) Velocity dispersion measurements for a sample of 374 Type 2 Seyfert galaxies from Heckman et al. (2004), with z < 0.03 and S/N > 10 in 
the CaT region. We have measured cr* using both the CaT and the Ca K regions. The solid line denotes (T(Ca K) = (T(CaT); the dashed line shows the fit de¬ 
rived from Figure 11 for an alternate set of Seyfert 2 galaxies (Eq. 2; see text for details), the dash-dot line shows the fit derived for this sample (Eq. 3), and 
the dotted lines mark the SDSS resolution limit. ([(T(Ca K) —(T(CaT)]/(T(CaT)} = 0.29 ±0.29 km s“’. (b) The same as (a), except that we have corrected the 
Ca K values using Equation 2, shown in (d) as a dashed line. We denote equality of the two velocity dispersion measures with a solid line. After the correction, 
{[<T(Ca K)-(T(CaT)]/t7(CaT)} = 0.11 ±0.34 km s-'. 


H(5 index tabulated by Brinchmann et al. (2004; Kauffmann et 
al. 2003a).‘* In the case of the larger sample the correlation, 
while highly significant, is weak (Kendall’s r = 0.22, Pnuii = 
10“^), while for the CaT comparison sample the correlation 
strength is even weaker (Kendall’s t = 0.14, Pnuii = 5 x 10““*). 
The largest values of (5 ct* correspondingly have the largest H(5 
indices, but a large H(5 index does not ensure a large . While 
young stellar populations apparently do play a role in driving 
the observed deviations, they may not fully account for it. We 
have also investigated whether including H(5 as a secondary pa¬ 
rameter improves the correlation between (T(Ca K) and (t( 4000- 
6200 A) or (T(CaT), using the partial Kendall’s r test from Akri- 
tas & Siebert (1996). In both cases the partial correlation coeffi¬ 
cient is lower than the correlation between the two cr* variables 
alone. We thus do not have a complete explanation for the offset 
at this time. 

Nevertheless, we can use the trend observed in Eigure 11 to 
derive a correction for the Ca K velocity dispersions. We use an 
ordinary least-squares fit, and derive the errors using bootstrap 
simulations. We find 

cr* = 1.40(±0.04)cr(CaH±K)-71(±5). (2) 

We also calculate a correction using the CaT comparison sam¬ 
ple only: 


The two corrections are not consistent within the estimated 
uncertainties, but when we apply Equation 2 to the CaT 
data set, we do significantly improve the agreement between 
cr(Ca K) and cr(CaT). Prior to correction, ((5 ct*) = ([cr(Ca K)- 
cr(CaT)]/cr(CaT)) = 0.29 ± 0.29. After applying the correction 
from Equation (2; see Eig. \2b), we find ((5 ct*) =0.11 ±0.34, 
while when we use the correction derived from this data set, 
we find ((5 ct*) = 0.0015 ±0.35. Clearly the correction derived 
from the CaT data is superior in the sense that the mean differ¬ 
ence is closer to zero after the correction is applied. However, it 
is worth noting that the scatter increases significantly (~ 20%) 
in both cases. Although the scatter increases, thus limiting the 
achievable precision, the most crucial point is that we can con¬ 
trol the systematic bias in the case of the Type 2 objects shown 
here. 

If this is a valid correction, we should be able to apply it to 
our AGNs, a completely independent sample of objects, and 
improve their agreement with CaT. We have performed this test 
in Eigure IQb, using Equation 2. In general we are overcorrect¬ 
ing at low cr*, leading to (Jcr*) = -0.15 ±0.39 (the scatter and 
offset become worse when Eq. 3 is employed). To determine 
whether Type 1 AGNs in general can be treated in this manner 
will require a larger database of velocity dispersion measure¬ 
ments spanning a wider range in cr*. We are currently assem¬ 
bling just such a sample (J. E. Greene & L. C. Ho, in prepa- 


cr* = 1.7(±0.1)cr(CaH±K)-121(±13). (3) 

The H<5 stellar absorption index gives a measure of the age of the stellar population; younger populations yield a larger H<5 index. 










Measuring cr* in AGNs 


15 



Fig. 13.— Input-output simulations of velocity dispersion measurements from the Ca H+K region, similar to Figure 5, but with higher AGN contamination. 
Galaxy models constructed as the average of 10 template stars. = [(j(out) —cr(in)]/cr(in). Simulations were carried out for S/N from 15 to 55, in steps of 5. The 
highest S/N (i.e., no additional noise) is 95 (shown in open circles). As in Figure 5, we highlight in solid lines S/N = 30 {triangles) and S/N = 50 {boxes). We 
show an AGN component comprising {left) 90% and {right) 98% of the total continuum. 


ration). While corrections of this nature are very dangerous, 
depending as they do on the detailed stellar populations of a 
given sample, we emphasize again that the Ca H+K lines may 
be the only available stellar features for estimating cr* in high- 
luminosity AGNs. To illustrate that the Ca K line can be used 
even at very high levels of AGN dilution, in Figure 13 we show 
simulations for AGN continuum fractions of 90% and 98%. As 
can be seen, even with an AGN fraction as high as 90% the 
Ca K line can still provide an estimate of cr* to within a factor 
of 2, for a moderate S/N of 20. However, with 98% AGN dilu¬ 
tion, there is no hope of obtaining a meaningful measurement 
of cr*. 

5. PRACTICAL GUIDELINES 

The optimal spectral region for measuring cr* depends on the 
Eddington ratio and continuum level of the AGN, as well as the 
redshift of interest. For z < 0.05, CaT is the region of choice. 
In the redshift interval of 0.05 < z < 0.76, the Fe region is the 
best option, provided that the Eddington ratio is ^0.5 and the 
overall dilution level is ^85%. Otherwise, the Ca K line must 
be used. Finally, for 0.76 < z < 1.3 Ca K is the only avail¬ 
able choice (assuming the observations are made in the optical 
band). 

We summarize these findings schematically in Figure 14. For 
a given Mbh, we estimate the bulge luminosity using the rela¬ 
tion of Marconi & Hunt (2003), log (Mbh/Mq) = 1.26(±0.13) 
[log (Lb/Lq)-10] -2.96(±0.48). If we make the simplifying 
assumption that the entire bulge and only the bulge is contribut¬ 
ing to the starlight, we can then estimate the AGN dilution in 
the B band for any AGN luminosity (expressed here in Edding¬ 
ton units). Of course, this is an approximation that is aperture- 
dependent. At z = 0.1, the projected SDSS aperture of 5 kpc 
is well matched to the effective radius /?eff ~ 4 kpc of a cr* = 
150 km s“' bulge (Bernard! et al. 2003; Greene & Ho 2005a), 
while at higher z we are more likely underestimating the to¬ 


tal contribution from starlight. We further caution that the ac¬ 
tual AGN contamination will differ between spectral regions, 
becoming worse toward the blue (see §4.1). Finally, there is 
significant intrinsic scatter in the MBH-Tbuige relation; Marconi 
& Hunt estimate that the intrinsic scatter is 0.48 dex. Using 
corrections from Fukugita et al. (1995) to convert from mg to 
m,, we further indicate the maximum redshift to which an AGN 
of m, = 19.1 mag (the SDSS magnitude limit for AGN-selected 
spectroscopic selection) would be detected by the SDSS. The 
lines of constant redshift shift to the right or left for brighter 
or fainter apparent magnitudes, respectively. Bear in mind that 
these lines are schematic only. We do not account for the com¬ 
bined (AGN plus host galaxy) color of the system, and it would 
be more appropriate to apply the SDSS galaxy magnitude limit 
below a certain AGN fraction. However, our goal here is simply 
to illustrate how the different constraints define natural regions 
in the Mbh-/agn plane. 

We have assumed that the observations are limited to opti¬ 
cal wavelengths. Although spectroscopy in the near-infrared is 
challenging, due to telluric absorption features and strong and 
variable emission from the night sky, it may be the best option 
in some cases. For instance, narrow-line Seyfert 1 galaxies, 
thought to be low-mass BHs at high Lboi/TEdd (e-g-, Pounds et 
al. 1995; Boiler et al. 1996), typically contain such strong Fe II 
contamination that CaT may be the only measurable absorption 
line. Of course, the chances of strong Ca II emission are also 
highest in these objects, making the determination of cr* partic¬ 
ularly challenging (e.g., Botte et al. 2004; Barth et al. 2005). 

6. SUMMARY 

We present a statistically significant sample of 40 AGNs with 
measurable velocity dispersion (cr*) using regions surrounding 
the Ca II triplet, the Mg lb triplet, and Ca H+K lines, for the 
purposes of intercomparison. Using a newly developed direct- 
fitting code, we perform a comprehensive set of simulations and 
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log / M J 

Fig. 14.— Schematic diagram to delineate the regimes in which the dif¬ 
ferent spectral regions can be most effectively used to measure cr*. The solid 
curves mark the three key redshift intervals (z = 0.05, 0.76, and 1.3) to which 
an mi = 19.1 mag AGN will be spectroscopically targeted by the SDSS, and 
the dashed curves identify, from bottom to top, Tboi/^Edd = 0.01, 0.1, 0.5, 
and 1. The fraction of light coming from the bulge was estimated using the 
MsH-f^buige relation of Marconi & Hunt (2003). 

cross-checks to evaluate the merits and limitations of using 
each spectral region, with the aim of obtaining realistic esti¬ 
mates of the many systematic uncertainties that affect measure¬ 
ments of (T* in AGNs. We argue that the Ca II triplet is least 
susceptible to template mismatch and AGN contamination from 
emission lines, and so provides the most reliable measurements 
of (j*. We therefore use these lines as a benchmark to test the 
other spectral regions. 

We examine two types of AGN contamination: featureless 
continuum dilution that effectively lowers the S/N of the ab¬ 
sorption features, and emission lines, both narrow and broad, 
that fill in the stellar absorption lines and bias the line profiles in 
subtle ways. In terms of dilution, we find that cr* is measurable 
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for AGN fractions ^71%, ^85%, and ^90% using Ca II triplet, 
Mg lb, and Ca Hh-K, respectively. As for AGN emission-line 
contamination, we find that measurements around Mg lb are 
very sensitive both to broad Ee II emission and narrow emis¬ 
sion from [Ee VI], [Ee VII], and [N I]. The narrow lines can be 
avoided by using the spectral region 5250-5820 A. Broad Ee II 
contamination can bias the fits severely when the Ee II width is 
narrowest. Eor a given Mbh this will occur at the highest lu¬ 
minosity, and thus is worst when the object radiates at a large 
fraction of its Eddington luminosity. In such cases, and for the 
highest levels of AGN dilution, the Ca Hh-K spectral region 
(and particularly the Ca K line) is the best remaining option. 
While there is a systematic offset between a(Ca K) and de¬ 
rived from other spectral regions, in Type 2 AGNs it is well fit 
by a linear relation that allows us to derive unbiased ct* esti¬ 
mates from Ca K. Eurther work is required to ensure that this 
correction is generally applicable. 

As a by-product of this study, we have derived reliable ve¬ 
locity dispersions for 35 Type 1 AGNs, and set out reasonable 
selection criteria to generate a much larger sample of AGNs 
with a broader range of z, black hole mass, and accretion rates. 
Such a sample would enable us to examine the Mbh “ ct* re¬ 
lation for active galaxies with true statistical power. We will 
be able to calibrate virial masses for black holes with unprece¬ 
dented accuracy, as well as search for second-order trends in 
the Mbh “ relation with mass, accretion rate, and redshift. 
Greene & Ho (2005c) present a preliminary investigation of the 
Mbh - f7* relation using the sample analyzed in this study. 
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